CONTROLLABILITY OF ROLLING WITHOUT TWISTING OR 
SLIPPING IN HIGHER DIMENSIONS 



ERLEND GRONG 



Abstract. We describe how the dynamical system of rolling two n-dimensional 
connected, oriented Riemannian manifolds M and M without twisting or slip- 
ping, can be lifted to a nonholonomic system of elements in the product of the 
oriented orthonormal frame bundles belonging to the manifolds. By consider- 
ing the lifted problem and using properties of the elements in the respective 
principal Ehresmann connections, we obtain sufficient conditions for the local 
controllability of the system in terms of the curvature tensors and the sec- 
tional curvatures of the manifolds involved. We also give some results for the 
particular cases when M and M are locally symmetric or complete. 

1. Introduction 

The rolling of two surfaces, without twisting or slipping, is a good illustration of 
a nonholonomic mechanical system, whose properties are intimately connected with 
geometry. It has therefore received much interest, and we can mention [T| 14} [5l ITT] 
and [21 Chapter 24] as examples of research produced in this area. In particular, 
the treatment of rolling in [21 [4] was done by formulating it as an intrinsic problem, 
independent of the imbedding of the surfaces into Euclidean space. 

The generalization of this concept to that of an n-dimensional manifold rolling 
without twisting or slipping on the n-dimensional Euclidean space, is well known 
(see e.g. [321 p. 268], [101 Chapter 2.1]). It is usually formulated intrinsically, in 
terms of frame bundles, and is an important tool in stochastic calculus on mani- 
folds. A definition for two arbitrary n-dimensional manifolds rolling on each other 
without twisting or slipping, first appeared in [201 App. B], however, this only dealt 
with manifolds imbedded into Euclidean space. An intrinsic definition for rolling of 
higher dimensional manifolds, that connected the definition in [20] with the intrin- 
sic approach in [21 S] , was presented in [8] . Apart from appearing as mechanical 
systems, rolling of higher dimensional manifolds can be also used as a tool in in- 
terpolation theory. For demonstration of the "rolling and wrapping" -technique, we 
refer to [11]. See also [21] for an example where this is applied in robot motion 
planning. 

For the rolling of two 2-dimensional manifolds, there is a beautiful correspon- 
dence between the degree of control and the geometry of the manifolds [21 E]. 
Essentially, we have complete control over our dynamical system if the respective 
Gaussian curvatures M and M do not coincide. Controllability in higher dimen- 
sions has been addressed in some special cases [TBI [22] ■ The first general result 
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on controllability in higher dimension is presented in [6], where it is shown that 
the curvature tensors determine the brackets of the distribution obtained from the 
constraints given by neither twisting, nor slipping. It is also proved that being able 
to find a rolling to an arbitrary close configuration, connecting the same two points, 
is a sufficient condition for complete controllability. 

The objective of this paper will be to describe the connection between geome- 
try and controllability for rolling of higher dimensional manifolds. We do this by 
connecting the earlier mentioned viewpoint from stochastic calculus with the one 
presented in [5]. 

This paper is organized as follows. In Section 2, we state the intrinsic defini- 
tion of rolling. We present some of the theory of frame bundles, and develop our 
notation there. We end this section by showing how we can lift our problem to 
the oriented orthonormal frame bundles of the involved manifolds. We continue in 
section Section 3, by doing computations on the lifted problem. By using properties 
of the sections in the Ehresmann connections, we obtain formulas for computation 
of the brackets of the rolling distribution. In Section 4 we project the results back 
to our configuration space of relative positions of the manifolds. Section 4.2 con- 
sist of conditions for controllability in terms of the Riemann curvature tensor and 
the sectional curvature of the manifolds involved. We end this section with some 
examples. Section 5 focuses on results concerning the rolling of locally symmetric 
and complete manifolds. Section 6 contains a brief comment on how to general- 
ize the concept of rolling without twisting or slipping to manifolds with an affine 
connection, and why the results presented here also holds for a rolling of manifolds 
with a torsion free affine connection. 

The author would like to express his gratitude to Mauricio Godoy Molina for 
many fruitful discussions concerning this subject. 

2. Intrinsic definition of rolling and its relations to frame bundles 

2.1. Intrinsic definition of rolling without twisting or slipping. Throughout 
this paper, M and M will denote connected, oriented, n-dimcnsional Ricmannian 
manifolds. Since in the special case n = 1, the conditions of rolling without twisting 
or slipping become holonomic (see [8]), we will always assume that n > 2. We adopt 
the convention to equip objects (points, projection, etc.) related to M with a hat 
("). Objects related to both of them are usually denoted by a bar ("), while objects 
connected to M are not given any special distinction. The exception to this rule is 
the Riemannian metric and the affine Levi-Civita connection which are respectively 
denoted by (•, •) and V on both M and M. The context will make it clear which 
manifold these objects are related to. For a vector field X on M, we will write X\ m 
rather than X(m), and we will use similar notation for other sections of bundles. 

For any pair of oriented inner product spaces V and V, we let SO(V, V) denote 
the space of all orientation preserving linear isometries from V to V. This allows 
us to define the SO(n)-fiber bundle Q over M x M by 

Q = |q e SO (r m M, TfoM\ : m e M, m e m\ . 

We can be sure that this fiber bundle is principal in the case when n = 2, but not 
in general. The space Q represents all configurations or relative positions of M 
and M, so that the two manifolds lie tangent to each other at some pair of points. 
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The isometry q : T m M — > Tf^M, represents a configuration where M at m lies 
tangent to M at m. The relative positioning of their tangent spaces is given by 
how q maps T m M into T^M . A rolling then becomes a curve in the space of these 
configurations. 

Definition 1. Let it and n denote the respective natural projections from Q to M 
and M . A rolling without twisting or slipping is an absolutely continuous curve 

q : [0,t] -> Q, with m(t) := n(q(t)), fh{t) := n(q(t)), 

satisfying the following conditions: 

No slip condition: m(t) — q(t)rh(t) for almost every t, 
No twist condition: an arbitrary vector field X(t) is parallel along m{t) 
if and only if q(t)X(t) is parallel along fh(t). 

From now on we will mostly refer to a rolling q(t) without twisting or slipping 
as simply a rolling. 

These two conditions can be described in terms of a distribution D of rank 
n on Q. Consider the problem of finding the rolling q(t) between two different 
configurations qo and q\ such that m(t) (and thereby also fh(t)) has minimal length. 
Here, by a distribution on Q, we mean a sub-bundle of TQ. This type of curves can 
be viewed as optimal curves in an input-linear drift free optimal control problem or 
length minimizers in a sub-Riemannian manifold. We will describe this structure 
in more detail in Section 12.31 First, however, we will review some facts about 
connections and frame bundles. 

2.2. From principal bundles to oriented orthonormal frame bundles. Con- 
sider a general principal G-bundle t : P — »• M, where the Lie group G acts on the 
right. We call the sub-bundle V := kerr* of TP the vertical space of P. If g is the 
Lie algebra of G, then for any element A £ g, we have a vector field a (A) defined 

by 

(1) a(A)\ p( f>= ± <Kp-exp G (tA)), foranypeP^eC^P). 
at t-o 

Here, exp G : g — >• G is the group exponential. We remark that a (A) is a section of 
V and for any p £ P, the map g — > V p , A n> a(A)\ p is a linear isomorphism. 

A sub-bundle £ of TP is called a principal Ehresmann connection if TP = 
£ © V and satisfy r g £ p = £ p . g , where r g denotes the right multiplication of g £ G. 
Equivalently, we can consider a principal connection form, which is a g-valued one- 
from satisfying the two conditions 

r*uj = Ad((? _1 )w, and u)(a(A)\ p ) = A for any A £ Q,p £ P. 

There is a one-to-one correspondence between these two structures, in the sense 
that kerw is a principal Ehresmann connection when u is a principal connection 
form, and for any principal Ehresmann connection £ , we can define a principal 
connection form by formula 

oj(v) = A for v £ T p P if v - a(A)\ p £ £ p . 

Related to a choice of Ehresmann connection £ we also have horizontal lifts, 
since the mapping r*|£ p : £ p — > T T ^P is a linear isometry. For a vector v £ T m M, 
we define the horizontal lift h p v of v at p £ r _1 (m) as the unique element in £ p 
which is projected to v by r*. 
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Let us consider a particular principal bundle over a manifold M. For two vector 
spaces V and V', let GL(V, V) be the space of all linear isomorphisms from V to 
V. For any m £ M, we say that a frame at m is a choice of basis /i, . . . , /„ for 
T m M. Equivalently, we can consider a frame as a linear map / G GL(R™,T m M). 
The correspondence between the two point of views, is given by 

/(0,...,0,1,0,...,0) = / J , 

V v ' 

1 in the jth coordinate 

We write F m (M) = GL(M n , T m M) for the space of frames in T m M. There is a 
natural action of GL(n) on these spaces by composition on the right. Using this 
action, we define the frame bundle f : J-{M) ->Mas the principal GL(n)-bundle 
with fiber over m being F m {M). 

Similarly, if M is an oriented Ricmannian manifold, we can define the oriented 
orthonormal frame bundle r : F(M) — > M as the principal SO(n)-bundle whose 
fiber is F m (M) := SO(M n , T m M). Here, R n is furnished with the standard orien- 
tation and the Eucliean metric. 

Let V be an affine connection on M, seen as an operator on vector fields 
(X, Y) i-> VxY- Then we can associate a principal Ehresmann connection £y 
on t : F(M) — > M to V, by defining £y to be the distribution on J~(M) consisting 
of tangent vectors of smooth curves f(t) such that the vector fields fi(t), . . . , f n (t) 
are all parallel along m(t) := f(f(t)). If M is Riemannian and oriented, and if V 
is compatible with a metric, then £y can be defined as a distribution on F(M) 
instead, since both positive orientation and orthonormality are preserved under 
parallel transport. 

We now go to the concrete case where V is the Levi-Civita connection. Define 
the tautological one-from 9 — (6i, . . . , 9 n ) on F(M), as the R"-valued one-form 

6i\t=T*Wjl t:F(M)^M, 

where b : TM — > T*M is the isomorphism induced by the Riemannian metric. 
In other words, if v G TfF(M), then Oj(v) = (fj,T*v). Denote the so(n)-valued 
principal connection form corresponding to the Levi-Civita connection by u>. The 
formulas for the differentials of 9 and lo are given by the well-known Cartan equa- 
tions. We express them in notations, that will be helpful for later purposes. 
Let R be the Riemann curvature tensor, defined by 

R(Y U Y 2 ,Y 3 ,Y 4 ) = (R(Y U Y 2 )Y 37 Y 4 ) , where R(Y U Y 2 ) = V^V^-V^V^-Vp^]. 

Furthermore, use it to define the curvature form £1 = (Oy), as the so(n)-valued 
two- form 

(2) ny(«i,«2) = R{T*V U T*V 2 ,fj,fi), V U V 2 G T f F(M). 

Then the following relations hold 

d °j + E"=i ^3* A °i = o, 

(3) 

dlOij + X/fe=l ^ ^kj = Qij, 

where uiij are the matrix entries of us. These equations are going to be important 
in order to understand the connections between geometry and the control system 
of rolling manifolds. 



CONTROLLABILITY OF ROLLING IN HIGHER DIM. 



5 



Remark 1. The Cartan equations can also be denned on the frame bundle corre- 
sponding to a general affine connection. See, e.g., [14 [ [20 ] for details. 

From now on, we will adopt the convention that whenever we mention K™, it will 
always come furnished with the standard orientation and the Euclidean metric. 

2.3. The rolling distribution and the corresponding sub-Riemannian struc- 
ture. The tangent vectors of all possible rollings form an n-dimensional distribu- 
tion D on Q. We will call this distribution D the rolling distribution. A curve q{t) 
is a rolling, if and only if, it is horizontal with respect to D, i.e. it is absolutely 
continuous and q(t) G A?(t) for almost any t. 

We present the following local description of D (see [8] for more details) . Let us 
write the projection of Q to M x M , as W : Q — > M x M. Given any sufficiently small 
neighborhood U on M , let e be a local section of the oriented orthogonal frame 
bundle F(M) with domain U . We write this local section as (e, U). Let (e, U) be a 
similar local section of F(M). We can use these local sections to trivialize the fiber 
bundle Q over U x U , with the map 

q E SO(T m M, T ffi M) i ^ (to, to, (q tj )) eU xU x SO(n), q t] := (fi^qe-). 

Then the distribution D\ w ^ 1 ^ Ux fjy in the above coordinates, is spanned by the 
vector fields 

(4) ej := ej + qe 3 ; + {(ea^e.ep) - (qe a ,\7 qej qep)) W^. 

l<a<P<n 

where j = 1, . . . ,n. Here, e 3 is seen as a vector field on Q\ Ux jj = 7f _1 (C7 x U) and 
qej stands for the vector field q H> q{ej\^^ q )). The vector fields W^g are defined by 




The symbol I here is not a parameter; it simply stands for "left" (an explanation 
of this will follow in Remark [5]) . 

Let us consider the optimal control problem finding of a rolling q{t) connecting 
two configurations qo and qi, such that the curve m(i) = Ti{q{t)) has minimal 
length. We do this by introducing a metric (•, •} on D, defined by 

(6) (V1,V2) = (7T*tJl,7IvU2), Vl,V2 € D q . 

With respect to this metric, the ej in (j4]) are a local orthonormal basis. Also, from 
the definition of D, we have (vi,V2) — (7r*tJi, 7r*W2}- Hence, minimizing the length 
of m(t) is equivalent to minimizing the length of fh(t) = Tr(q(t)) (this is can also be 
seen from the no-slip condition). 

Definition 2. A triple (Q, D,{-, •)), where Q is a connected manifold, D is a 
distribution on Q and (■,■) is a metric on D, is called a sub-Riemannian manifold. 

The sub-Riemannian distance function d(qo, qi) between two points is defined as 
the infimum of the length of all curves which are horizontal to D. We can view 
rollings between configuration from go to q\ along a curve of minimal length as a 
sub-Ricmannian length minimizer in Q. 

For the distance function to be finite, we need that every pair of configurations 
qo and q\ can be connected by a curve horizontal to D. In other words, we need to 
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determine when rolling without twisting or slipping is a controllable system. Our 
goal is to give sufficient conditions for this to hold, in terms of geometric invariants 
on M and M. We will leave the question of finding optimal curves for later research. 

Remark 2. The vector fields W*p in (J5J) can be considered as a "locally left invari- 
ant" basis of ker7r». It will be practical to also introduce a "locally right invariant" 
analogue. Relative to two chosen local sections (e, U) and (e, U), define 

(7) ^=g^ s _l__ gQS _^_), qir .= {e uqej ). 
Notice that W r afi = Ya, s Qaiqi3 s Wf s . 

2.4. Controllability and brackets. Given an initial configuration qo £ Q, write 
O qo for all points in Q that are reachable by a rolling starting from go- This will 
be the orbit of D at qo, which coincides with the reachable set of D, since D is 
a distribution (see, e.g., [3J [TS] for details). The Orbit Theorem [UJ [52] tells us 
that O qo is a connected, immersed submanifold of Q, but also that the size can be 
approximated by the brackets of D. Define the C°°((3)-module Lie D as the limit 
of the process 

D 1 =T(D), D k+1 = D k + [D, D k ], 
T(D) denotes the sections of D. Let D k and Lie g D be the subspaces of T q Q 
obtained by evaluating respectively D k and LieD at q. Then, for any q G O qo , 

(8) Lie 9 D C Tq O qo . 

In particular, it follows from ([5]), that if D is bracket generating at q, i.e., if 
Lie 9 D = T q Q, then O qo is an open submanifold of Q, and we say that we have 
local controllability at qo- If O qo = Q for one (and hence all) go G Q, the system is 
called completely controllable. 

The least amount of control happens when O qo is rt-dimensional submanifold. 
As a consequence of the Orbit theorem and Frobenius theorem, this happens if and 
only if D\q is involutive, that is, if Lie g D — D q , for every q £ O qa . 

The focus of this paper will be to provide results of controllability, by investi- 
gating when the distribution D will be bracket generating at a given point q. 

Remark 3. When D is not bracket generating at q, Lie q D will in general only give 
us a lower bound for the size of O q . However, if LieD is locally finitely generated 
as a C 00 (Q)-module, i.e., has a finite basis of vector field when restricted to a 
sufficiently small neighborhood, then the equality holds in ([5]). 

Remark 4. We will use the notation introduced here for distributions in general, 
not just for the rolling distribution. 

2.5. Relationship between frame bundles and rolling. Consider the linear 
Lie algebra so(n). For integers a and /? between 1 and n and not equal, let w a p be 
the matrix with 1 at entry a/3, -1 at entry /3a, and zero at all other entries. Clearly 
w a p = —wp a . Define w aa — 0. The commutator bracket between these matrices 
are given by 

The collection of all w a /3 with a < /3 form a basis for so(n). 
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Consider the principal bundle F(M) — > M. We introduce the following notation. 
Write <7q/3 for the vector fields on F(M) corresponding to w a p in the sense of |T]). 
Let £ be the principal Ehresmann connection corresponding to the afhne Levi-Civita 
connection on M, with corresponding principal connection form oj. As before, we 
denote the tautological one-form by 9. By using horizontal lifts with respect to £ , 
we can define vector fields Xj on F{M) by 

(10) Xjlf-hffj, 

which satisfy 9i(Xj) = Si.j. Hence, the tangent bundle of F(M) is trivial, since it is 
spanned by {Xj}" =1 and {cr a/ 3} Q <^. Finally, we define d a p , X 3 ■, Q and 6 similarly 
on F(M). 

The configuration space Q may be identified with F(M) x F(M) quotiented out 
by the diagonal action of SO(n). Let w denote the principal SO(n)-bundle 

m : F(M) x F(M) ->• F(M) x F(M)/SO(n) = Q. 

Then w(f, f) = q, if / = q o f. By viewing uj, D, 9 and 9 as forms on F{M) x F{M), 
we are able to obtain the following result. 

Theorem 1. Let 

V = kerw n kerw n ker(9 - 9), 

and let D be the rolling distribution. Then m*T> = D, and the map is a linear 
isomorphism on every fiber. 

Proof. From its definition, it is clear that {Xj + is a basis for V. 

Choose any pair of local section (e, U) and (e, U) of respectively F(M) and 
F(M). Give F(M) x F(M)\ U q local coordinates by associating the pair of frames 
(/, /) to the element 

(11) (m,m, (/,,•), (£,•)) eUxUx SO(n) x SO(n), 

if fj = Z)"=i fij e i\m and /y = Yh=i fijZi\m holds. 

Relative to this trivialization, we can define left and right vector field on each of 
the SO(n)-factors. On the first, define 

(12) ^ = 5( /sQ ^" /s ^)' K " = ^{ f0s ^~ fas ^)' 

Notice that ^l r a a = Y^i s =i falf ps^u- Remark also that ^ l a p is just the restriction 
of cr Q( 3 to F{M) x F(M)\ Ux jj, while ^ r a g depends on the chosen local section e. 
Define "3/^ and W a p analogously on the second SO(n)-factor. 

Restricted to F{M) x F(M)\ Ux q and in the local coordinates (TTT1) . the vector 
fields Xj and Xj can be written as 

n n 

(is) Xj = e u (e. - E . ^ = E & («» - E f *A) . 

s=l a<,9 i-l a</3 

where = (e Q , V £i e^) and = (e a , Ve 4 e». 
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We now turn to the image of T(F(M) xF(M)\ Ux q) under zu*. Define qtj , ej, 
and W r a p on Q\ Ux Q as in Section 12751 Remark [2] allows us to rewrite ej on the form 



e, = e s + qe 3 + £ (Tf^Wi - £ q sj T^W^ 

a</3 8=1 

Locally the mapping w can be described as 

n 

■ (m, fh, (fij), (fijfj h-> (m, fh, (%.,•)) , q {j = f is fj. 



w 



and the action on the tangent vectors is given by formulas 

w * : * r fl h-> -W l R " 

ap f * r ap 

$r i — ^ W r „ 

ap ' r r ap 

From this it is clear that V = D and that the map is injective of each fiber, 
since 



= uj* 22 fjs \ Xs + X, 



8=1 



□ 



From the form of the distribution T>, we obtain the following interpretation of 
rolling without twisting or slipping. 

Corollary 1. Let q(t) be a rolling without twisting or slipping. Let (f(t),f(t)) be 
any lifting of q(t) to a curve in F(M) x F(M) that is horizontal to T>, and define 
m(t) and fh(t) as the respective projections to M and M. Then {f{t),f(t)) satisfy 
the following 

(i) (No twist condition) Every vector field fj(t) is parallel along m(t). Every 
vector field fj(t) is parallel along m(t). 

(ii) (No slip condition) For almost every t, 

f-\t)(fa{t))=r\t)(fh{t)). 

Furthermore, if (f(t), f(t)) is any absolutely continuous curve in F(M) x F(M), 
satisfying (i) and (ii), then w(f(t),f(t)) is a rolling without twisting or slipping. 

Proof, (i) follows from the definition of the principal Ehresmann connections ker u> 
and kerw. (ii) is exactly the requirement for a curve to be in ker(# — 0). □ 

The main advantage of the viewpoint given in Theorem [TJ is that it will help us 
to compute Lie 9 -D. 

Corollary 2. ro* T>\ = D k q for any q G Q, k 6 N. 

Proof. We only need to show this locally. Introduce local coordinates as in the 
proof of Theorem [1] and let ej be as in Q . Then 

n n 

J2fis (x s +X s ),J2frs (; V - ' V 

,8=1 8=1 
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and since X)"=i fis (x a + X s ^j is a local basis for £>, it follows that D 2 = vu* T> 2 
locally. The rest follows by induction. □ 

Since [JQ, Xj] — 0, computations of brackets of T>, and hence of brackets D, can 
be reduced to mostly computing brackets of sections in the Ehresmann connections 
corresponding to the Levi-Civita connections of the manifolds involved. 

2.6. Remark on previous descriptions of rolling using frame bundles. The 

description of rolling given in Theorem [TJ looks very similar to the definition of 
rolling without twisting or slipping found in [4] for dimension 2. Here, the descrip- 
tion of a rolling is in terms of the distribution T> := £>0ker?u*, which can also be 
described as 

(14) V = kcr(w - w) n ker(6» - 6). 

In [4], a rolling of a pair of 2-dimensional manifold is defined as a curve in Q, that 
is horizontal to zu^V, where V is defined in terms of (1141). 



We could have used V for our computation, since clearly zu*V q = D q also, and 
D is bracket generating at a point q £ Q if and only if T> is bracket generating at 
any (and hence every) (/, /) G vj~ 1 (q). However, since [V k , kertu*] c V k + ker m* , 
the additional brackets are not of any interest. 

The definition of rolling or "rolling without slipping" in probability theory is 
defined on frame bundles [10j . and is equivalent to considering curves in F(M) x 
F(M) that are horizontal to T> for the special case when M is K™ with the Euclidean 
metric and standard orientation. 

3. Brackets of V 

3.1. Tensors on M and associated vector fields. We introduce a general type 
of vector fields associated to tensors. By giving a general formula for the brackets 
of these, we essentially determine all the brackets of V. Let T be a tensor on M. 
We will only consider tensors defined on vectors. To any tensor fc-tensor on M, we 
can associate the functions 

g iu ... ih (T): F(M) -> R 

/ ^ T(f n ,...J lk ) ' ^*.^»- 

If k = 2 + 1, and T is antisymmetric in the first two arguments, we can define vector 
fields on F(M) by 

l<a</3<n 

If T is a tensor on M, we define S'i 1 ,...i h (T) and similarly as respec- 

tively functions and vector fields on F(M). 

Lemma 1. Let Xk be defined as in (|10ll . 

(a) For any I -tensor T , 

Xk£ii,...,u(T) = <£ S ii,...,i ! ,fc(VT). 

(b) For any 2 + I -tensor T , that is antisymmetric in the first two arguments 

n 

[X k , W iu ... til (T)} - W iu ..., iuk (VT) - J2 ^.,fc,u,..,<. (T)Xs- 

s=l 
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Proof. (a) We use a local representation of Xj . Consider the formula given in 
(|13p . and write X)T=i fsj e s & s just fj- Since ^ a p is just the representation 
of o~ a p in local coordinates, we can write Xj as 

l<a<n 

In this notation, remark first that o~ a pfi — $p,ifa — &a,ifp-> which gives us 

1 - 

X (fa, Vf k fp)<Tapfij =2 X! Uai'J f k ffi) {5f3,ija ~ 5 a ,ijp) 
l<a<P<n a,P=l 

n 

=E(^> v /^>^ = v /*4- 

Then the result follows from realizing that 
Xk&i u ...,i l (T) = fkTifii, . . . - 2J (fcf,'Vf h fp)o- a pT(f il ,...,fi l ) 

l<a<P<n 

= fkT(fi t , • - • , /i, ) - T( V/ fc /<! , . . . , /j, ) - TX/ij , V/ fc /is, • . . , /i, ) 

T(/ il ,/ i2 ,...,V /(c /i i ) 

= VT(/i, , . . . , /ij,/fc). 

(b) The brackets [er Q ^, ct^a] are, by definition, given by the same relations as 
described in ([9]). We will continue by the following computations. 



k , A— 1 a,(3,K, A— 1 

^ n 1 n 

— 2 XI ^K,X,H,-,il( T ) a ^fk + J X/ ^«,A,ii,...,i ( (7 1 )o-«;A ((/a, ^f k fp))^aP 
k,A— 1 a,/3,K,A— 1 

i " 

I - 

+ 2 E ( T (^f k f^fx,fi 1 ,---,fi l ) + T(f K ,Vf h fx,fi 1 ,...,f il ))a K> , 

K,A=1 

n ^ n 

-E^,fe,i 1 ,...,i ! (T)/ s + - ^ &8, a ,i u ...,i,(T)((f a ,Vf k fp})<T a p 

s=l a, 0,8=1 

I " 

+ 2 E ^.,fc,n,-,«(T)«/a.V/.^»«7 QiS 

a,/S,s=l 
1 - 

+ 2 X ^s,/3,U,...^( T ) V>- i A ((/ Q , V/ fe / 5 .}) Cr Q/3 



a,^,s=l 

'ii,...,ii,fc(VT) - y^^ > a ,fc,i 1 ,...,ii(T)-X'a 
s=l 



□ 
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The next lemma gives an explanation for the introduction of the above notation. 
Lemma 2. [X h Xj] = - W l2 (R) for i, j = 1 . . . n. 

Proof. This lemma is an easy consequence of the Cartan equations. Since ker# n 
kerw only contains the zero section of TF(M), we can show equality in the above 
equation by evaluating the left and right hand side by 9 and oj and see that it 
produces the same result. Evaluating the left hand side, we get 

e{[x u x 3 ]) = -de{x i ,x j ) = o, 

w([Xi, X,]) = -duiXuXj) = J2 ^(X^Xj) 

l<a<f3<n 

l<a<0<n 

which is obviously what we get from evaluating the right hand side. □ 

Remark 5. Combining these two lemmas, we get a way to express the commuta- 
tors of the Ehresmann connection. Let £ = kerw be the Ehresmann connection 
corresponding to the Levi-Civita connection. Then 

gk+2 = ^^{^^fl)}^ , 

for k > 0. As a consequence of this, we get the well know Ambrose-Singer theorem, 
see [3],[lH App C], that the subalgebra spanned by elements w|/ (/^ / i 1 ,...,i fc+2 (V fe i?)) 
is contained in the homology algebra at /. 

We adopt the convention that if the elements in the collection are vector fields, 
"span" means the span over smooth functions (so in Remark it means over 
C°°(F(M))). If the elements are vectors, the span is over the real numbers. 

3.2. Obtaining the brackets for T). Computing the brackets on V is a bit more 
complicated than each individual Ehresmann connection, since it is harder to know 
whether or not two vectors fields are equal modspan{Xj + Aj}™ =1 rather than just 
modspan{X,-}™ =1 . We illustrate this by computing the two next brackets. 

Lemma 3. (a) [X k , [X U X 3 ]] = - W ijk {VR) + £Li S skij (R)X S . 
(b) Let R 2 be the 6-tensor on M defined by 

R 2 (Y a , Y P , Y n , Y l2 , Y l3 ,Y U ) = R(R(Y a , Y p )Y il , Y l2 , Y l3 ,Y U ). 

Then 

[X u [X k , [Xi, Xj]] - - W l]k i{V 2 R) + W lkl] {R 2 ) 

n 

H" ^ 1 ijslk 

{VR)-£ ijskl (VR))X s . 
The reason for the notation R 2 will be clearer in Section [5Tl 
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Proof. Statement (a) follows directly from Lemma [TJ By Lemma [T] we also have 
[XuiXk^X^Xj]]] 

n n n 

= - #^ ; (V 2 i?) + *Mik(VR)X a + J2 <? S kiji(VR)X a + J2 * a uj(R) W U R ) 

S—l S—l S — l 

n 

= - Wijki{V 2 R) + W mj {R 2 ) + fe*fc( Vi *) - <?ijski(VR))x s . 

8=1 

□ 

We can continue this procedure, computing more of the brackets using Lemma[T] 
However, these will become more and more complicated. Also, for a general pair of 
manifolds, it is hard to determine which brackets actually give us something new, 
that is, something that could not be expressed as linear combinations of previously 
obtained vectors. Rather than giving the total picture, we will therefore focus on 
giving some sufficient conditions, which are usually more simple to check. 

4. Sufficient conditions for controllability 

Let R and R be the curvature tensor on M and M respectively. Define a new 
4-tensor of elements in D, by 

R = TT*(R) -Tt*{R). 

Remark that R may also be seen as a bilinear map of two elements in f\ 2 D. Use 
Vi? to denote the 5-tensor on D, defined by 7r*(Vi?) — 7f*(VR). Finally, introduce 
a bundle morphism 1Z : /\ D — > /\ D* , so that 

2 

(15) ^(!i)(? 2 )=S(? 2) ? 2 ) ^ 2 e/\D q . 

4.1. Projection of the results on V. From the discussion in previous section, 
we have the following formulations for the brackets of D. 

Lemma 4. Let (e, U) and (e, U) be two local sections of F{M) and F(M), respec- 
tively. Then on Q\ Ux fj, in terms of the notation introduced in Section \2.!A 

D 2 = D 1 span j ^ ep, e h e^W^X 

^ l<a<f3<n ' *>J =1 



D 3 = D 2 + spani ^ VR(e a ,e l 3,e i ,e j ,e k )W^ l3 + qR(ei,ej)e k - R(qei,qej)qe k 

l<a<0<n 

• U k=l 

Proof. The formula for D 2 follows directly from Lemma[5]and the local formulation 
of nj, given in the proof of Theorem [T] To see how the expression for D 3 follows 
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from Lemma [3] (a) , observe first that 

n n 

^ [£ s kij (R)X S + £ s kij (R)X s ^j = — ^ ($ ijks (R)X S + S ijks (R)X f 

s=l s=l 
n 

= {fa* (• R ) - # « *» (A))** ( mod v ) ■ 

3=1 

Furthermore 

n 

fipfjxfkK (^fiXKs(R) — <^hXks(R) \ X s 

s,/i, A,k, — 1 

= qR(ei,ej)e k - R(qei,qej)qe k - ^ (qe a , V qR{e ^ ej)ek qep) W l afi 

l<a</3<n 

n 

l<a<P<n 



□ 



Corollary 3. Define a bundle morphism S : D © /\ 2 D — > £>* fey 

2 

=%oft(£), veD q ^e/\D q , 

where b- : rj ^ rj(v, •) /or any zT 6 Z? 9 and two /orm 77. TTien dimD^ > n + 
rank f2 1 g + rank S 1 9 . 

Proof. Given point g g Q, introduce local coordinates in a neighborhood of q, in 
the way demonstrated in Section 12731 We will also keep the same notation from the 
previous mentioned section. Then, all we need to show is that for a given q G Q, 
the dimension of 

(16) span{gi?(e i ,e i )e fc (7r(g)) - R{qei,qej)qe k (Tf(q))}^ j>k=1 C T^ q) M, 
is equal to rankS 9 . 

Introduce isomorphisms t> : D — >• £>* and jj : £>* — > £>, relative to the metric 
defined in Use the same symbols for the isomorphisms between the tangent 
bundle and the cotangent bundle on M and M, induced by their respective metrics. 
Observe that 

S(e fc ,e, A e,j) = tt* {^{R{e i ,e j )e k )^ -tt* (\>{R{qe u qej)qek\ 
= \>Y - \>Z, 

where 

Y = R(ei,ej)e k + qRfa, e 3 )e k 

+ (( e «' ^ R{eu ei )e k ep) - (qe a , V qR ( eue j) ek qep)) W^, 

l<a<P<n 

Z =q~ 1 R(qe i ,qe J )qe k + R(qei,qej)qe k 

+ (( e «' V g - 1 %e i , g e J )ge fe e ^}-(9 e «' V % ei ,ge,-) g e fc 9 e /?)) W t&- 

l<a</3<n 
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From this, it becomes clear that 7r*j} is a bijective linear map from the image of S q 



4.2. Sufficient condition in terms of the curvature tensor and sectional 
curvature. As mentioned before, there is a strong connection between controlla- 
bility and geometry in the two dimensional case. 

Theorem 2 ( 2, 4]). For q £ Q, let x q denote the Gaussian curvature of M at 
n(q), and let 5c q denote the Gaussian curvature of M at n(q). Then 



If k — >t = on O q , then dim O q — 2. 

The "if and only if in the above theorem follows from the fact that in two 
dimensions, the rolling distribution D at a point q, is either bracket generating or 
involutive. This does not hold in higher dimensions, however, but we are able to 
present the following generalization. 

Definition 3. The smallest integer k such that = Lie 9 D is called the step of 
D at q. 

Theorem 3. Let 1Z be as defined in (|15l) . Then, for any element q,q$ S Q, the 
following holds. 

(a) dimO go — n if an only iflZ\o m = 0. 

(b) IfTZ q is an isomorphism, then D is bracket generating of step 3 at q. 

Hence, if O qo contains a point q, so that lZ\ q is an isomorphism, then 
O qo is an open submanifold. 

Remark 6. The statement in Theorem[3] (a) was also presented in [6l Cor. 5.28]. By 
combining 6, Cor. 5.26] with [51 Prop. 5.17], and doing some simple calculations, 
we can also obtain the result of Theorem 3(b), however, this is not stated. The 
proof is presented here, since the approach in [B] differ from ours, and since the 
results were obtained independently. 

Proof. Statement (a) becomes obvious from Lemma |4j To prove (b), let 7f(g) = 
(m,m), where If : Q — > M x M is the projection. Pick respective local sections 
(e, U) and (e, U) of F(M) and F(M) around m and m and use them to introduce 
the local coordinates defined in Section [2. 31 Since lZ\ q is an isomorphism, we know 
that 



to (HU). 



□ 



dim O q 



= 5, 



if and only if 



x — >c ^ on O q . 





Lemma |4] then tells us that 



D\ = span{ ej (q),qei(m), W r ^(g)}" J> ^=i 

= spanjej^.e^m), W^(g)}^ )Q;/ 3 =1 = T q Q. 



□ 
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To state that lZ q is an isomorphism is equivalent to claiming that R q induces a 
pseudo-inner product on /\ 2 D q , i.e. a nondegenerate bilinear map. Therefore, we 
have the following way we can check if 1Z is an isomorphism at q. 

Corollary 4. Choose any orthonormal basis {vj} of T m M,m = ir(q). Compute 
the determinant of the "("^^ x "C^" 1 ) matrix 

det (R(v a ,vp,Vi,Vj) - R(qv a ,qvp,qvi,qvjf), 

1 < a < (3 < n are row indices, 1 < i < j < n are column indices, 
If this determinant is nonzero, then we have local controllability at q. 

From this we obtain the following corollary. 

Corollary 5. Define a function Jc q on 2- dimensional planes L in D q by the formula 

X q {L) = X^ {q) {lT^L) - X^ {q) (TT,L), 

where >c m and Stf^ denotes the respective sectional curvatures of M and M at the 
indicated points. Then 

(a) dim O qo = n if an only if ~x q = for any q G O qo . 

(b) IfJc q >0 or~K q < 0, then D is bracket generating of step 3 at q. 

Proof. li~>c q = 0, then R q is also. Similarly, if H p > (resp. >? p < 0) for every L, 
then — R q (resp. R q ) will be an inner product on /\ 2 D q . 

To see this, for the case Jc q > 0, we only need to show that R q (£, £) < 0, whenever 
£ G /\ 2 D q is nonzero. Pick an orthonormal basis basis vi, . . . ,v n of D q . Write 

Lij := span{vi,Vj}. 

In this basis, we have that if £ = J2i<i<j< n a ijVi A Uj is nonzero, then 

l<i<j<n 

The case Jc q < is treated similarly. □ 

Remark 7. All the conditions stated here, are sufficient conditions for local con- 
trollability. However, if they hold in all points, they will naturally be sufficient 
conditions for complete controllability. 

4.3. Examples. 

Example 1. We start with two known examples, to verify our results and demon- 
strate their effectiveness of obtaining information on controllability. 

(a) If M is a sphere of radius r and M — R" is the n dimensional Euclidean 
space, then M has constant sectional curvature \, while M has constant 
sectional curvature 0. It follows that ~x q = \ for any q E Q. Hence D is 
bracket generating at all points, and the system is completely controllable. 

(b) If M and M are the spheres with respective radii r and f, then 

_ _ 1 1 

for any q G Q. Hence the system is completely controllable if and only if 
r 7^ r . When r — r, D becomes an involutive distribution. 



16 



E. GRONG 



To compare, see [HI [23] for a former proof of the controllability of (a), and [16] for 
a treatment of the example in (b). 

Example 2. More generally, if M is any manifold with only strictly positive or 
strictly negative sectional curvature, rolling on n-dimensional Euclidean space, then 
this system is completely controllable (we will later show that this only needs to 
hold in one point of M). 

Example 3. Let M = S 2 x 1 be the subset the Euclidian space K 4 , 

{(xo,xi,X2,Xi) el 4 : xl + x\ + x\ = 1} . 

Define a local section on the subset U = {(xo,Xi ) X2,Xs) G M : x 2 > 0}, with the 
orthonormal vector fields 



ei = -Jx\ + xl ( -d xo + 2 X ° 2 {xid Xl +x 2 d X2 ) 
V x i + X 2 

e 2 = , * 2 : ( ~d xi + —d X2 ) , e 3 = d X3 . 
y/x{ + xj V X2 J 

(a) Let us first consider M rolling on R 4 . The rolling distribution can locally 
be describes by 

D 1 = span{ei, e 2 , e 3 }. 

Xf) 

ei=ei+gei, e 2 = e 2 + qe 2 H -===W 12 , e 3 = e 3 + qe 3 . 

V X 1 + X 2 

D is then of step 3 for any q G U and 

D 2 = D 1 8 span{Wi 2 }, D 3 = D 2 © spanjgei, qe 2 } 

Since D 3 is locally finitely generated, we know that dim O q = 6 for any 
q G U (and for symmetry reasons, every q G Q). 

(b) Let M roll on a copy of itself. Consider the rotation matrix (qij) = 
((e{, qej}). Give (q^) the coordinates 

(cos 9 cos y> sin 9 cos i/> — cos 9 sin sin i/> sin 9 sin i/j + cos 9 sin cos ip 
— sin cos ip cos cos i/ 1 + sin 9 sin y cos ip cos sin ip — sin # sin ip cos t/j 
— sin ip — sin i/j cos tp cos <p cos */? 

The vector fields spanning D are locally given by 

_ xo (sin cos "0 + cos sin <^ sin -0) 
d = ei + qei V, 

V Zi + ^2 

- x o tta a^o (cos 6 cos + sin sin w sin 

e 2 = e 2 + ge 2 + VKi 2 V, 

y xf + £ 2 V ^l + x 2 

_ xo cos w sin ip 
e 3 = e 2 + qe 2 H V, 

V x i + X 2 

V := cos ip cos V'W / i 2 — cos ^ sin ^>Wi 3 — sin ipW 23 . 
The matrix — fi?(e a , ep, ~e~i, ej)j, i < j,a < ft is then given by 

1 — cos 2 ip cos 2 ip — cos 2 if sin */j cos t/j cos ip sin y cos ip 
— cos <p sin ip cos ^ — cos ip sin 3/1 cos sin ip sin i/ 7 

cos ip sin <£> cos i/ 1 cos ip sin <p sin ip — sin 2 <£> 
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It is easily checked that this matrix has rank 2, except when sin ip = sin ip = 
0. Restricted to the subset of Q where the latter equation hold, that is, the 
configurations where the two copies of the line 

A = {(0,0,0,x 3 )GM}, 

lie tangent to each other, D is involutiove and the orbits are 3 dimensional. 
On the other points, we have that 

D 2 = D 1 ©span{Wi2, F}, D 3 = D 2 © span{qei, qe 2 , qe 3 }. 

so the orbits have dimension 8, or codimcnsion 1. 

This example illustrates that if we are rolling manifolds of dimension 
higher than two, the dimension of O q does not only depend on the connect- 
ing pair of points Tf(q) = (m,fh). 

5. Particular cases 

We present some special results for when the manifolds involved in the rolling 
are particular nice. We will first deal with locally symmetric spaces, then present 
some results for rolling on complete spaces. Remark that all of these results are 
applicable to the case of rolling on R™, since this is both locally symmetric and 
complete. 

5.1. Locally symmetric spaces. Recall the definition of f2 from ([2]). 

Proposition 1. Let M be locally symmetric and let M be flat (R = 0). Then D 
is at most of step 3. D is bracket generating at q G Q if and only if 

2 

is a linear isomorphism. 

Proof. Consider the bundle V, and let Xj and Xj be defined as in (fT0|) . Since M 
is flat [Xi + Xi, Xi + Xj] = [Xi,Xj]. Then 

n 

[X l7 X 3 ] - - Wij{R), [X k , [X u X,]] = £ ski j{R)X s , 

s-1 

n 

[X h [X k , [Xi,Xj]]] =-^2^ skij {R) W sl (R) g V 2 , 

s=l 

Hence, V 2 f ^ + kerw <f ^ = T {f ^ (f(M) x F(M)) only if 

span{r y (i2)| (/)f) }". =i = B Pa n{a i j\ (f>?) Y i ^ i . 

and this also implies that spanfe^, 2 s =i ^skij{R)X s } = spanjcrij, fj\, which gives 
us the desired result. □ 

When M is not flat, the results become a little bit more complicated, and require 
us to introduce some notation. Let R l denote the 21 + 2-tensor defined by R 1 = R 
and 

R l (Y ai Yp,Y h ,Y i2 ,..., Y i2l _ 1} Y 2l ) := R^WY*, Yp)Y h ,Y i2 ,..., Y hl _, , Y 2l ). 
Lemma 5. // VR = 0, then WR l = for any 1>1. 



18 



E. GRONG 



Proof. We give the proof by induction. Assume that VR k = for 1 < k < I. Let 
m(t) be any smooth curve in M. Let vx(t), . . . , v n (t) be parallel vector fields along 
mit). Then 

^i? J (fQ!,W/3)Uu,«i3,-< -,V2l-l,V2l) 

d n 

= -j7 2J R(Va,Vf3,V n , V s )R l ~ 1 (v s ,V i2 V 2 l-l,V 2 l) = 0. 

s— 1 

Hence VR l = also. □ 

We introduce a notation related to R l , similar to what we did for R. Use R l 
for the analogues tensor on M, and write R for the tensor on D defined by R = 
TT*(R l )-n*(R l ). 

Proposition 2. Let M and M both be locally symmetric. Then D is bracket 
generating at q if and only if 

(17) |Jspan J2 ^( 

e ai ep, ejj , . . . , ei 2l )W a /3 \ q \ ~ ker 7r* \ q . 

1>1 l<a</3<n > 

Proof. We will look at the brackets of V. From Lemma El we know that for any 
I > 1, 

[X h , [X i2 , [■■■ [X i2l _ 1 ,X l2l ] ■■■]= (-1)' ^h,...,i 2 ,( R )> 

n 

[x n , [x l2 , [■ ■ ■ [x l2l ,x i2[+1 ] • • •] = (-i) /+1 2J $ s,ii,-.j2i+i (R l )x s - 

s=l 

Analogues relations hold for the brackets of Xj. Projecting the even brackets to 
T q Q, we get the left hand side of (IT71) . which has to be equal to all of ker7f*| g in 
order for D to be bracket generating. Conversely, if (fT7|) holds, then the projection 
of the odd brackets will span T q Q together with ker7f*| g and D q . □ 

5.2. Rolling on a complete manifold. The fact that one of the manifolds is 
complete, makes it easier to give statements about complete controllability. The 
reason for this, can be summed up in the following Lemma. 

Lemma 6. Assume that M is complete. Let t 1— ¥ m(t) be any absolutely continuous 
curve in M with domain [0,r]. Let qo € Q be any point with Tt(qo) — m(0). 
Then there is a rolling 1 1— > q(t) of M on M, so that 

q{0) = qo, no q(t) = m(t) for any t 6 [0, r]. 

Proof. If we assume first that both M and M are complete, then such a rolling q(t) 
exist. The proof for this can be found in [2j p. 386]. This proof is done for the 
case when M and M are two dimensional, but can, with simple modifications, be 
generalized to higher dimensions. 

Assume now that M is not complete. Let f(t) be a lifting of m(t) to a curve 
in F(M) that is horizontal to the Ehresmann connection, i.e. each fj(t) is parallel 
along m(i). Define the curve in M™ by, 

m(t) = [ f- 1 {s)(m(s))ds. 
Jo 
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Let f(t) be a lifting of m(t) to a curve F(M. n ), so that each fj(t) is parallel along 
fh(t). Then, from CorollaryCO q{t) = w(f(t), /(*)) is a rolling of M on E™. 

Let go := 9(0) • Since both R" and M are complete, we know that there is a 
rolling q(t) of M™ on M along m(t), so that g(0) = go So" 1 - We can then obtain 
our desired rolling by defining q(t) = q(t) o g(t). □ 

Proposition 3. Let the manifold M be complete. Assume that there is a point 
m G M, so that D q is bracket generating for every point q G 7T (m). Then the 
system is completely controllable. 

Proof. Let go be any element in Q. From Lemma [5] we know that there is a rolling 
q(t) from go to some point gi G n~ l (m). Hence O qa = O qi . But since D is bracket 
generating in gi , O qi is an open submanifold. Since go was arbitrary, we have local 
controllability at every point, so O qo = Q for any go G Q. □ 

Corollary 6. Let M be a manifold that is both complete and flat. Assume that 
there is a point m G M, so that for some ( and hence any) orthonormal basis {v^J—i 
ofT m M, 

det(R(v a ,vp,Vi,Vj)) ^ 0. 
1 < a < ft < 1 are row indices, 1 < i < j ' < 1 are column indices. 
Then the system is completely controllable. 

Example 4. Let M be the surface in K 3 , defined by 



M = {(x 1 ,x 2 ,x 3 ) GM 3 : ^xj+xj =l-f(x 1 ),\ Xl \ < §}, 

where 



if|xi|<l 
L e-Clxxi-i)" if 1 < j^l < 3 
Define the following orthonormal basis on M, 



C2 = 77 — z{~x 3 d X2 +x 2 d X3 ). 

1 - f{xi) 

All Christoffel symbols are determined by 



T{ 2 = (ei, V ei e 2 ) = 0, Tl 2 = (ei, V e2 e 2 ) = 



(W(x 1 )) v /l + /'(^i) 2 
and from this we can compute the Gaussian curvature 

_ f'jx,) 

^ (l + /'(x 1 )2)2(l_/( a;i ))- 

Inserting the value of fix) we obtain that k[x) = for \x\\ < 1, but strictly positive 
for 1< \x x \ < §. 

It follows that , if we roll M on R 2 , the system is completely controllable. Observe 
that in this case 

Lie£> = spanjei + gei,e 2 + ge 2 + T\ 2 Wi 2l f{x\)W 12 , f{xi)qe u f{xi)qe 2 } , 
fails to be locally finitely generated around points fulfilling \x\ \ = 1. 
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6. Further generalization of rolling without twisting or slipping 

Up until now, we have only been concerned with rolling two Riemannian mani- 
folds on each other without twisting or slipping. The definition can easily be gen- 
eralized to manifolds with an affine connection. We introduce the generalization 
here. 

Let M and M be two connected manifolds, with respective affine connections V 
and V. Then a rolling without twisting or slipping can be seen as an absolutely 
continuous curve q(t) into the manifold 

Q=[qE GL(T m , T m M) : m G M, m G m\ . 

satisfying (no slip condition) and (no twist condition) from section \2. II 

Reexamining the proofs, it turns out that the description of many of the results 

we had for rolling related to the Levi-Civita connection, generalizes to general 

connections. We will describe this here briefly. 

Define the tautological one-form on J-(M) as the R n -valued one-form obtained 

by 

0{v) = r 1 {v), veT f F{M). 

Define 9 similarly, while U) and Q are defined in terms of the connection. We will 
also choose Riemannian structures on M and M in order to present the results 
in a similar way and make them easier to compare, but these do not need to be 
compatible with the connections. 

The rolling distribution D, will still be an n dimensional distribution, and the 
relation in Theorem [1] is still valid. D is locally spanned by vector fields 

n 

ej = ej + qej + ^ ((e Q , W ej ep) - (qe a ,V 'qe^e^jj E\ r 

a,P=l 

where =E s =i9«gf-. 

The study of controllability becomes harder, since torsion may appear in the 
equations ([3]). However, this is actually the only complication. Let e a p G jjl(n,R) 
be the matrix with only 1 in at entry a/3 and zero otherwise. Define Xj as in 
equation flu]) . If we then modify the definition of 7P'i 1 ,...,i l (T) to 

n 

"^i u ...,it{T)= £<*,p t h,...,i,(T)o-(e a p), 

a, 0=1 

now defined for any tensor T, then Lemma [T] still holds for any connection. Since 
all of our results follow from Lemma [1] and [2J it follows that our results also holds 
for any pair of manifolds with torsion free connections. 
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